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Vacancy complexes with oversized impurities in Si and Ge 
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In this paper we examine the electronic and geometrical structure of impurity- vacancy complexes 
in Si and Ge. Already Watkins suggested that in Si the pairing of Sn with the vacancy produces a 
complex with the Sn-atom at the bond center and the vacancy split into two half vacancies on the 
neighboring sites. Within the framework of density-functional theory we use two complementary ab 
initio methods, the pseudopotential plane wave (PPW) method and the all-electron Kohn-Korringa- 
Rostoker (KKR) method, to investigate the structure of vacancy complexes with 11 different sp- 
impurities. For the case of Sn in Si, we confirm the split configuration and obtain good agreement 
with EPR data of Watkins. In general we find that all impurities of the 5sp and 6sp series in Si 
and Ge prefer the split-vacancy configuration, with an energy gain of 0.5 to 1 eV compared to the 
substitutional complex. On the other hand, impurities of the 3sp and 4sp series form a (slightly 
distorted) substitutional complex. Al impurities show an exception from this rule, forming a split 
complex in Si and a strongly distorted substitutional complex in Ge. We find a strong correlation 
of these data with the size of the isolated impurities, being defined via the lattice relaxations of the 
nearest neighbors. 

PACS numbers: 71.55.Cn,76.60.-k,61.72.-y,61.72.Ji 
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I. INTRODUCTION 

Intrinsic defects and their complexes with impurities 
play an important role in semiconductor physics. In 
this paper we present ab initio calculations for vacancy- 
impurity complexes in Si and Ge. Based on EPR mea- 
surements, Watking-'^^ has already shown that in Si the 
Sn-vacancy complex prefers an exotic configuration, i.e. 
a split-vacancy complex with the Sn atom on the bond- 
center position and the vacancy split into two half- 
vacancies on the nearest neighbor (NN) sites. This con- 
figuration has been supported by ab initio calculations 
for Sn in Si by Nylandsted Larsen et al^ and Kauko- 
nen et ali^, reporting a small energy preference of 0.045 
eV for this complex compared to the configuration with 
substitutional Sn and the vacancy on NN sites. 

Recently we have studied by density functional calcu- 
lations the Cd-vacancy and Cd-interstitial complexes in 
Si and Ge^, aiming at understanding the electric field 
gradients (EFG) measured in perturbed angular corre- 
lation experiments for the ^^^In/^^^Cd probe atom. We 
find that both in Si and Ge the substitutional Cd-vacancy 
complex is instable and relaxes into the highly symmetri- 
cal Cd-split-vacancy complex, being about leV lower in 
energy than the substitutional configuration. In the split 
configuration the Cd atom hybridizes very weakly with 
the six Si or Ge nearest neighbors, resulting in a nearly 
isotropic charge density and an extremely small EFG. 
In this way we find good agreement with PAC measure- 
ments and can uniquely assign the very small measured 
Cd EEC's of 28 MHz in Si and 54 MHz in Ge to the 
Cd-split-vacancy complex. In parallel to our calculations 
Alippi et. al. show in a recent publication'^, that in Si 
also In impurities form such a split configuration with the 
vacancy, exhibiting an unusually large binding energy of 
2.4 eV. 



The present paper has a twofold aim. First we present 
ab initio calculations for the electronic and geometrical 
structure of the Sn-vacancy complex in Si and Ge. In 
particular we concentrate on the hyperfine properties of 
this defect for the different charge states, i.e. the hyper- 
fine fields, the isomer shifts and the electric field gradi- 
ents, and compare with available experimental data by 
Watkinsi and othersSii. These properties have not been 
calculated so far, since they are difficult to obtain by the 
standard pseudopotential plane wave method. In the sec- 
ond part of the paper the structure of impurity- vacancy 
complexes for other impurities is calculated, in particular 
heavy impurities with larger sizes than Si or Ge atoms. 
For the elements Cd, In, Sn and Sb of the 5sp-series we 
find that both in Si and Ge the split-vacancy impurity 
complex is preferred over the substitutional one by ener- 
gies between 0.5 to 1 eV, being slightly larger in Si than 
in Ge. The results for the even heavier element Bi of 
the 6sp series suggest, that in Si and Ge this is the sta- 
ble configuration for all oversized impurities. In order to 
discuss the importance of impurity size for the relative 
stability of the substitutional and split configurations we 
calculate the lattice relaxations of the nearest neighbors 
for a large series of isolated impurities from the 3sp, 4sp, 
and 5sp series. In general, we find a good qualitative cor- 
relation of the sign and size of the NN relaxations with 
the stability of the two configurations, although a strict 
one-to-one correspondence is not valid. 



II. THEORETICAL METHODS 

All calculations are based on density functional the- 
ory in the local density approximation^. Two different 
methods have been used to solve the Kohn-Sham equa- 
tions. The first one is the pseudopotential plane wave 



(PPW) method, which has mostly been used to inves- 
tigate the configuration and stabiUty of the impurity- 
vacancy configurations. These configurations, as well as 
the relaxations of the neighboring atoms, have been recal- 
culated by the KKR-Grccn- function method^' ^° for the 
Sn-impurity, which as an all-electron method allows also 
to calculate electric field gradients, isomer shifts and hy- 
perfine fields^ ^■^^. 

The PPW method approximates the inhomogeneous 
systems containing defect complexes by periodically ar- 
ranged large supercells, and uses plane waves to ex- 
pand the electronic wave functions. This has the ad- 
vantage that band-structure methods can be used to 
determine the electronic structure, and total energies 
and forces on the atoms can be calculated without dif- 
ficulty for arbitrary arrangements of the atoms in the 
supercell. We have used norm-conserving Kleinman- 
Bylander (KB)-pseudopotential&^^ for all atoms consid- 
ered, where the s,p- valence electrons are treated by pro- 
jectors, and the d{l — 2)-component is used as a local po- 
tential. Our EStCoMPP-program— is fully parallelized and 
can efficiently handle supercells with up to 300 atoms. 
For the calculations of impurity-vacancy complexes we 
used a (lll)-oriented supercell with the basis vectors 

61 = 3a(0,-l,l);62 = 3a(-l, 1, 0); 63 = 2a(l,l,l) con- 
taining 108 atoms, a is the theoretical lattice constant ( 
lO.Tlas for Ge with Ecut = 11.56 Ry, and 10.25aB for 
Si with Ecut = 9 Ry; Bohr radius qb = 0.529177 A). 
The impurity-complexes were placed in the middle of the 
cell. We used C31, symmetry explicitly for all configura- 
tions. The isolated substitutional impurities were calcu- 
lated with a 2x2x2 a"^ cubic cell containing 64 atoms. 
The impurities were located at the central site, and the 
Trf site-symmetry is enforced. A plane-wave basis set 
equivalent to 6 x 6 x 6 Monkhorst-Pack- A;-points was 
used which yields three inequivalent fc-points in the ir- 
reducible part of the Brillouin-zone for the 108-atom- 
supercell and five inequivalent fc-points for the 64-atom- 
supercell. We used a plane-wave cut-off Ecut — 9 Ry for 
impurities in Si, and Ecut = 11-56 Ry for impurities in 
Ge. The atoms belonging to the impurity-complexes and 
their nearest neighbors were relaxed until the forces on 
all atoms were less than 0.1 mRy/as. We checked that 
forces on further neighbors (which were not moved) were 
less than 20 mRy/as. The Sn-complexes in both hosts 
were recalculated for a cut-off energy of 13.67 Ry, and 
we found no significant changes of the configurations or 
the energy differences for the tested configurations. We 
estimate that the positions of all relaxed atoms are de- 
termined with an accuracy better than 0.01 a^, and the 
energy differences are accurate to about 0.1 eV. 

In the KKR-Green-function method the calculation is 
divided in two steps. First the Green function of the host 
is determined. In a second step the host Green function is 
used to determine the Green function of the crystal with 
a single impurity by a Dyson equation. For details we 
refer to referencei^. All calculations are performed with 
an angular momentum cut-off of Imax = 4. For the host 



crystal we use the LDA lattice constantsii of Si (10.21ao) 
and Ge (10.53ao). A k-mesh of 30 x 30 x 30 fc-points 
in the full Brillouin zone is used. We use the screened 
KKR-formalismi2ii2iiS with the TB-structure constants 
determined from a cluster of 65 repulsive potentials of 
4 Ry heights. The diamond structure is described by a 
unit cell with 4 basis sites, two for host atoms and two 
for vacant sites. For the Green function of the defec- 
tive systems, we allow 77 potentials of the defect and 
the surrounding host atoms to be perturbed, which are 
then calculated selfconsistently with proper embedding 
into the host crystal. All calculations include the fully 
anisotropic potentials in each cell and thus allow the reli- 
able calculation of forces, lattice relaxations and electric 
field gradients. The Coulomb potential can be expressed 
in terms of n{r) and therefore the force i^" on atom n can 
be derived with the use of the "ionic" Hellmann-Feynman 
theoremi^. The hyperfine parameters of interest in the 
present paper are the isomer shifts (IS), the Fermi con- 
tact term of the hyperfine field (HF) and the quadrupole 
splitting (A*). The last quantity is determined by the 
electric field at the probe-atom site. In a non-relativistic 
treatment the isomer shift can be calculated from the 
charge density n(0) at the nucleus 



IS ^ a n(0), 



(1) 



where a is the calibration constant. The hyperfine field 
is given by the Fermi contact contribution in terms of the 
magnetization density m(0) at the nucleus. 

AHF = ^fiB ™(0), (2) 

Finally the quadrupole splitting is calculated from 

A. = IfM , , ,3) 

^ ^source 

where Q is the quadrupole moment of the probe-nucleus 
and Vzz is the electric field gradient along the main (111) 
axis. The parameter c is the speed of light and Egource is 
the energy of the emitted 7-ray, which is 23.8 keV for Sn. 
The tensor of the EFG is given by the second derivatives 
of the Coulomb potential and can be written asi^: 



m=-2 
V2,miO) = h + h, 

5 io r'3 



(4) 



Rmt 



Rmt 



where l2m('^) are spherical harmonics, V2,mir) are the 
I = 2 components of the Coulomb potential and Rmt 
is the muffin tin radius of the probe atom. Furthermore 
?^2m(?') are the I = 2 components (m = —2, . . . , -1-2) of 




Figure 1: The Sn-split-vacancy configuration (left) and the 
substitutional Sn- vacancy configuration (right). 



the radial charge density. We note that just the / — 2 
components are needed for the EFGs, while the / = 1 
components enter into the forces. Consequently, the EFG 
vanishes for sites with cubic or tetrahedral symmetry. 
For the cell division of the crystal we used a generalized 
Voronoi construction. The impurity cell at the bond cen- 
ter was constructed such that it is slightly larger than the 
cells of the host atoms, thus avoiding to decrease the muf- 
fin tin radius of real atoms by more than 16%. For the 
treatment of occupied states in the band gap, we have 
used a separate energy contour for gap energies. Since 
we used group theory we were able to occupy these gap- 
states according to their symmetry, in this way avoiding 
to calculate the wave functions of the localized gap states. 



III. ELECTRONIC STRUCTURE OF 
Sn- VACANCY COMPLEXES 

As predicted already many years ago by Watkinsi on 
the basis of EPR measurements, and as has also be con- 
firmed by ah initio calculations?'?-, the stable configura- 
tion of the Sn-vacancy pair in Si consists of a Sn-atom 
at the bond-center position and the vacancy split in two 
halves on the two (empty) nearest neighbor positions. 
This split-vacancy configuration is shown schematically 
in Fig. ^together with the "normal" substitutional com- 
plex. As we have found recently^ the same split- vacancy 
configuration occurs also upon pairing with Cd-atoms, 
yielding both in Si and Ge very small electric field gra- 
dients, which allow to identify this complex uniquely. In 
this configuration the Sn or the Cd atom are coordinated 
by six host atoms at the relatively large distance of 1.26 
times the nearest neighbor distance (in the unrelaxed 
configuration). 

The local density of states (LDOS) of this complex at 
the Sn-site in the Ge host is shown in Fig. [3 (upper 
panel), together with the corresponding LDOS for the 
Cd split- vacancy complex in Go in Fig. |21 (lower panel). 
In the calculation the Fermi level is fixed at the maxi- 
mum E^ai of the valence band. The different peaks are 
labeled by the irreducible subspaces Axg, Ai^-, E^ and Eg 
of the Z?3rf group, representing the point symmetry of the 
complex. Compared to the Cd complex with the d-level 



at about -9 eV, the d-level of Sn is fully localized and at 
much lower energies. However otherwise the level scheme 
is very similar: the Aig state at relative low energies, the 
A2u and the doubly degenerate E^ state slightly below 
the Fermi level, and the doubly degenerate Eg state in 
the gap and, in the case of Sn, a second Aig state at 
higher energies. Since the occupied Aig, A2u and E^ 
states accommodate 8 electrons, of which the 6 neigh- 
boring host atoms provide 6 (one dangling bond each), 
the vacancy complex with Cd (with 2 valence electrons) 
is neutral for Ep — Eyai ([CdV]°), and the complex with 
Sn (with 4 valence electrons) is doubly positively charged 
([SnV]^+). As discussed in reference^, the level sequence 
is basically determined by the divacancy. The Cd and Sn- 
atoms can be considered as ionized Cd^^ and Sn^+ ions 
inserted in the center of the divacancy and only weakly 
hybridizing with the 6 nearest neighbors. The attractive 
ionic potentials shift the divacancy states to lower ener- 
gies, in particular the fully symmetrical Aig state, being 
the only state affected in the first Born approximation 
with respect to the potential. This effect is naturally 
stronger for the Sn''^+ ion than for the Cd^"*" ion, but the 
LDOS are very similar. Of course, also Cd and Sn s and 
p states, localized at energies above E^ah are hybridized 
into the occupied states, so that the local charge in the 
impurity cell is about 1 electron in the case of Cd in Ge 
and 2 electrons in the case of Sn in Ge. In fact, the Aig 
state above Ep can be considered as the genuine Sn 5s 
states, which is, however, not occupied. 

By occupying the Eg state in the gap, we obtain 
complexes with additional electron charges, i.e. [SnV]"^, 
[SnV]0, ... and analogously [CdV]", [CdV]^-, ... com- 
plexes. Both the [SnV]+ and the [SnV]° complexes are 
magnetic with total moments of 1 and 2 fiB, and with 
small local moments of 0.045 fiB (0.008 /is) and 0.091 /is 
(0.016 /is) on the Sn site in Si (Ge). Although we have 
calculated only the three charge states [SnV]^+, [SnV]"*" 
and [SnV]°, calculations as well as experiments show that 
also the negatively charged [SnV]~ and [SnV]^~ states ex- 
ist. Since for the complex [SnV]° the Eg majority state 
is completely filled having a moment of 2 fXB , we expect 
that also the [SnV]~ state is magnetic with a total mo- 
ment of 1 fiB and local moments similar to the [SnV]^" 
complex, while the [SnV]^~ is non-magnetic. The exis- 
tence of these five charge states arises from the fact, that 
the Eg state is very extended and has moreover a nodal 
plane at the Sn-site so that Coulomb effects are very 
small. This also explains the very small local moments. 
The situation is quite similar for Cd complexes, except 
that the local moments are somewhat smaller in Si, e.g. 
0.024 fiB for [CdV]~. The same level sequences we ex- 
pect also for other 5sp impurities to occur. For instance 
for In, the charge state should be [InV]+, if Ep is fixed at 
Eyai, and the [InV]°, [InV]~ and [InV]^~ configurations 
should be magnetic. In all calculations the complexes in 
Si are very similar to the above results for the Ge hosts. 

The calculations yield sizeable relaxations of the neigh- 
boring host atoms towards the Sn atoms, which are prac- 
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Figure 2: LDOS projected on the irreducible subspaces of the 
Did group at the Sn site (above) and at the Cd site (below) 
for the impurity split-vacancy complex in Ge. The nearest 
neighbors are fixed on ideal lattice positions. 



tically radially symmetrical and increase with the charge 
state. For [SnV]2+ the NN Ge atoms relax by 2.9% of 
the NN distance, for [SnV]+ by 4.2% and [SnV]0 by 5.4%. 
In Si the corresponding relaxations are somewhat larger: 
6%, 6.2% and 6.4%. The relaxations arise from the at- 
traction of the neighboring electrons to the Sn'*+ ion and 
increase with occupation of the Eg gap states, since these 
states are dangling-bond like and mostly localized at the 
host atoms. Therefore the neutral [SnV]° state has larger 
relaxations than the [SnV]^+, which is opposite to the 
normal expectations. 



IV. HYPERFINE PROPERTIES OF 
Sn- VACANCY COMPLEXES 

Already Watkin^ has measured by EPR hyperfine 
fields at the neighboring Si atoms adjacent to a SnV com- 
plex, which he identified as the neutral [SnV]° complex 
with total spin S = 1. Table ^ shows the calculated 
hyperfine fields at the Sn atom and the six Si or Ge 
neighbors adjacent to the magnetic [SnV]+ and [SnV]° 
split complexes and a comparison with the experimen- 
tal valuesiiSa. First we note, that the field for [SnV]° is 
about twice the field of [SnY]"*". This is a consequence 
of the large extent of the Eg state, representing locally 
a small perturbation, so that the forces increase linearly 
with the occupancy. Second, for [SnV]'^ we obtain good 



Host Si Ge 

Hhf [kG] Hhf [kG] 

[V|Sn[V]+ 139 6.00 
[V|Sn|V]" 211 14.59 
EXP 241 (386MHz) 



Host 



Si 
Hhf [kG] 



Ge 
Hhf [kG] 



[V|Sn|V]+ -39.59 -44.46 
[VlSnlV]" -82.96 -94.19 
EXP -91.02 (-77.00MHz) 

Table I: The calculated hyperfine field (HF) at the Sn atom 
(above) and the nearest neighbors (below) for the Sn-split- 
vacancy in Si and Ge is shown. A comparison to the experi- 
mental value of Watkins- is given. 



agreement with the experimental values. Here we note, 
that our calculations are non-relativistic, so that the rel- 
ativistic enhancementii of the hyperfine field is missing, 
which would presumably improve the agreement. For 
the [SnV]~ complex we expect a similar hyperfine field 
as for the [SnV]+ complex, since the total moment is 
the same. We have chosen in our calculations the spin- 
polarization such, that the total moment has a negative 
sign (spin down occupation). As a consequence the lo- 
cal moments for all sites are negative, since the Eg state 
determines the magnetic properties alone. That means 
that mioc/ Hhf < at the Sn-site and vniocl Hhf > at 
the Si nearest neighbor site. The change of sign of the 
hyperfine field from Si to Sn can be explained due to the 
different exchange interaction mechanism between the d- 
like "spin down" Eg electrons and the "spin down" and 
"spin up" s-core-electrons at the Sn-site. In the first case 
the interaction is attractive, leading to a smaller s charge 
density at the nucleus, whereas in the latter case the sit- 
uation is opposite. In total the hyperfine field changes 
its sign in comparison to the Si site. In addition to the 
hyperfine fields we have calculated the isomer shifts and 
quadrupole splittings. These observables can be mea- 
sured by the Mossbauer spectroscopy. Measurements us- 
ing the ^^^Sn isotope in Si and Ge were performed by 
Weyer et. aliS^ in 1980, which however do not agree with 
our calculated results. From various private discussions 
we conclude, that in these early experiments clustering 
of Sn-impurities might have occurred, leading to more 
complicated structures. On the other hand, our results 
are in good agreement with unpublished measurements 
of R. Sielemann (private communication). 

In TableElthe results for Si and in Table llTTl the results 
for Ge are presented. The calibration constant a was de- 
termined by calculating n(0) for several Sn-defects, for 
which well known isomer shifts exist. Svaneet. ali^itSS al- 
ready calculated in a relativistic treatment an a of 0.092 
Oq mms~^ , which is in good agreement to ours with tak- 
ing the Shirley factor™ for Sn of 2.48 into account. A 



system 


n(f=0) 


IS EFG A* 




K7^] 


[mms^^] [MHz] [mms^^] 


a-Sn 


87759.82 


ref 0.00 0.00 


subst. Sn 


87758.81 


-0.2200 0.00 0.00 


[V|Sn|V]^+ 


87759.71 


-0.0257 -14.90 0.3875 


[V|Sn|V]+ 


87759.36 


-0.1076 -10.59 0.2713 


[V|Sn|V]0 


87758.99 


-0.1942 -5.62 0.1473 



Table II: Isomer shifts and quadrupole splitting for the Sn- 
vacancy complexes in Si. (a = 0.2178 Oq mms~^, Esource — 
23.8 keV, Q = -0.124 barn) 



system 


n(f=0) 


IS EFG A* 




K'] 


[mms-i] [MHz] [mms-^] 


Q-Sn 


87759.82 


ref 0.00 0.00 


subst. Sn 


87759.26 


-0.1400 0.00 0.00 


[V|Sn|V]^+ 


87760.61 


0.17 -12.05 0.31 


[V|Sn|V]+ 


87760.55 


0.16 -8.08 0.1938 


[V|SnlV]0 


87760.46 


0.14 -3.97 0.1027 



Table III: Isomer shifts and quadrupole splitting for the Sn- 
vacancy complexes in Ge. (a — 0.2178 qq mms~^, Esource ~ 
23.8 keV, Q = -0.124 barn) 
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Figure 3: Relaxations of the Sn- vacancy complex in Ge. Pre- 
sented are PPW calculations with a cut-off energy of 13.67 
Ry, where only the Sn-atom is relaxed (full lines) and in ad- 
dition the six NN of the complex are relaxed (dashed lines). 
On the left, the total energy is shown relative to the energy 
of the fully relaxed configuration with Sn at the bond cen- 
ter site. On the right, the corresponding Hellmann-Feynman 
forces on the Sn-Atom are shown. 



closer view on the isomer shifts for the Sn-spht- vacancy 
configuration (in relation to the Mossbauer Hne of a-Sn) 
shows a decreasing charge density at the nucleus with the 
occupation of the Eg state in the band gap. This behav- 
ior can be explained due to the screening impact of the 
d-like Sg-electrons at the Sn-site on the s-electrons re- 
sponsible for the IS. As in the Cd case {Qcd = 0.83 ) the 
Sn-EFG's are rather small, since hybridization with the 
neighbors is weak. It is obvious, that the EFG increases 
with occupying the Eg state, since Px and py charges give 
positive contributions to the EFG at the Sn-site. 



V. SUBSTITUTIONAL VERSUS 

SPLIT- VACANCY COMPLEXES WITH 

IMPURITIES 

In this section we describe the results for a whole se- 
ries of calculations for impurity-vacancy complexes in Si 
and Ge. We would like to find out for each of the 11 
investigated impurities, which complex, substitutional or 
split-vacancy, is most stable, and to draw a border line 
between these two complex families as a function of the 
"size" of the impurities. To achieve this aim, we have 
performed PPW calculations for vacancy complexes with 
the impurities Al, Si, and P of the 3sp series, with Ga, 
Ge, As, and Se of the 4sp series, with Cd, In, Sn, and Sb 
of the 5sp series and finally with 6sp element Bi. In all 
cases only the neutral state has been calculated. First 
we discuss in detail the results for the Sn impurity in Ge. 

Fig. 121 shows the total energy and the force on a Sn 
atom, if this atom is moved adiabatically, i.e. by full re- 



laxation of all other atoms, from the substitutional po- 
sition (indicated by "relative distance 1", measured in 
units of the NN distance from the vacant site) to the 
bond center ("relative distance 0.5"). The two curves re- 
fer to the cases where only the Sn-impurity is allowed to 
relax (full lines) and where in addition to the Sn atom 
also the six host neighbors are relaxed (dashed lines). 
While by relaxation of only the Sn atom an intermediate 
position of about 0.75 is obtained, in the case of full re- 
laxations the Sn impurity moves without energy barrier 
into the stable bond center position at 0.5. 

The results for all considered impurities in both Si and 
Ge are listed in Table Hvl The second column gives the 
force on the impurity in the unrelaxed substitutional po- 
sition, i.e. when the impurity and all neighbors are fixed 
at ideal lattice sites. The third column gives the energy 
difference between both configurations, such that for neg- 
ative values the split configuration is stable and for posi- 
tive values the substitutional one. The last column gives 
the exact position of the impurity as measured in units 
of the NN distance from the vacant site. 

The results show that both in Si and Ge the 5sp im- 
purities Cd, In, Sn, and Sb as well as the 6sp impurity 
Bi on the substitutional site experience a strong force, 
which pushes them into the "free space" at the bond cen- 
ter, thus gaining typically an energy of 0.5-1.0 eV. Note 
that the energy gains are very substantial, much larger 
than the value 0.045 eV reported in^. One would intu- 
itively classify these heavy impurities as oversized in Si 
and Ge. In line with this argument is the fact that the 
forces as well as the resulting energy gains are somewhat 
larger in Si than in Ge, due to the smaller lattice con- 
stant of Si. All other impurities, with the exception of 



Impurity 


force A£; = 


- Esplit — Esubst 


stable pos. 


in Si 


[mRy] 


[eV] 


in (111) 


Al 


-22.4 


-0.15 


0.50 


P 


-11.5 


-H.27 


0.95 


Ga 


-K5.37 


+0.22 


1.04 


Ge 


-21.7 


+0.30 


0.87 


As 


-41.7 


+0.82 


0.90 


Se 


-44.9 


+0.55 


0.92 


Cd* 


-85.2* 


-1.04* 


0.50* 


In 


-77.2 


-0.69 


0.50 


Sn 


-105. 


-0.83 


0.50 


Sb 


-129. 


-0.68 


0.50 


Bi 


-182. 


-0.92 


0.50 


Impurity 


force A_E = 


- Esplit — Esubst 


stable pos. 


in Ge 


[mRy] 


[eV] 


in (111) 


Al 


-0.96 


+0.03 


0.82 


Si 


-4.5 


+0.30 


0.88 


P 


-11.7 


+0.81 


0.93 


Ga 


-H14.6 


+1.11 


1.04 


As 


-44.1 


+0.41 


0.87 


Se 


-37.62 


+0.35 


0.91 


Cd* 


-69.6* 


-1.01* 


0.50* 


In 


-43.2 


-0.47 


0.50 


Sn 


-79.0 


-0.60 


0.50 


Sb 


-109. 


-0.64 


0.50 


Bi 


-152. 


-0.86 


0.50 



Table IV: Relaxations of impurity- vacancy complexes in Si 
{Ecut = 9 Ry, above) and Ge {Ecut = 11.56 Ry, below): The 
first column gives the forces on the impurity atom in the sub- 
stitutional positions, a negative sign means that forces are 
directed to towards the the bond center. The second column 
gives the energy differences between the fully relaxed con- 
figurations with the impurity at the bond-center and at the 
substitutional site. The last column gives the final position 
of the impurity in (111) direction; 0.5 is the bond center and 
1.0 is the substitutional lattice position. (* was obtained by 
the PAW method with Ecut = 20.25 Ry). 



Al in Si, prefer the substitutional configuration, however 
partly with large relaxations from the ideal position. For 
instance, the Ge impurity in Si is shifted by 13% of the 
NN distance in the direction of the vacancy and the As 
atom in Ge by a similar amount. Of all impurities, only 
the Ga impurity in both Si and Ge is shifted towards its 
remaining three neighbors, i.e. away from the vacancy. 
Therefore it is tempting to classify the Ga impurity as 
"undersized" . The behavior of the Al impurity is most 
surprising. As a member of the 3sp series it should be 
smaller than the isoclectronic Ga impurity of the 4sp se- 
ries. However the calculation seems to indicate a larger 
size than Ga. The force at the ideal lattice site drives 
the impurity in the direction of the vacancy, resulting in 
a split configuration for the Al-vacancy complex in Si, 
and in a strongly distorted substitutional configuration 
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Figure 4: Nearest-neighbor relaxations for various 3sp, 4sp, 
5sp and 6sp impurities (substitutional point defects) in Si 
(above) and Ge (below). 



in Ge. 

Unfortunately, we do not know any reliable rules for 
the "size" of impurities. Pauling'^^ has given some general 
rules about the volume of atoms, e.g., in a tetrahedral en- 
vironment. According to this, the size of impurities of the 
same row of the periodic table should decrease with in- 
creasing valence, However, for the investigated cases we 
find just the opposite, i.e. the size increases with increas- 
ing valence, at least for the impurities of 4sp and 5sp 
series. For instance, As is larger than Ga, and Sb larger 
than In. Moreover, according to Pauling^i, the size of 
isovalent impurities should increase with increasing main 
quantum number. As expected, our calculations confirm 
this second rule in general. However, deviations from 
both rules occur for Al: We find that Al is larger than Si 
and P, and in particular, is larger than the isovalent Ga 
impurity, which belongs to the higher 4sp series. 

In order to shed some light on the size problem and its 
relation to the stability of impurity- vacancy complexes, 
we have calculated the lattice relaxations of the relevant 
isolated impurities in Si and Ge. For simplicity, we have 
only calculated the displacements of the nearest neigh- 



bors, by fixing all other atoms at the ideal positions. 
The idea behind these calculations is that both problems, 
i.e. the NN relaxations of the impurities and the different 
configurations of the impurity- vacancy complexes, are di- 
rectly correlated. For instance, the impurity exerts forces 
on the nearest neighbors, which, as a consequence, are 
displaced from their ideal positions. The impurity it- 
self is not shifted, since the reaction forces of the nearest 
neighbors cancel each other due to the tetrahedral sym- 
metry. However, in the case of a vacancy, one of the four 
neighbors is missing, so that the reaction forces or the re- 
maining three neighbors no longer cancel each other and 
shift the impurity either to or away from the vacancy. In 
this way the pair geometry should be related to the size 
of the single impurities. 

The calculated lattice relaxations of the nearest neigh- 
bors of the various impurities are shown in Fig. ^ They 
are given in percentages of the NN distance, with posi- 
tive values denoting outwards relaxations. Note that in 
the calculations only the NN are allowed to relax (the 
relaxation of the higher neighbor shells is expected to 
increase the NN relaxations, presumable by about 30 to 
50%). For the 4sp and 5sp impurities we clearly find that 
the relaxations increase with increasing valence of the im- 
purities. The same trend has been found in Table Hvl for 
the force on the unrelaxed impurity. Also the relaxations 
of the 5sp impurities are in general larger than the ones 
of the 4sp impurities, which is line with the fact that the 
forces on the unrelaxed impurity in the substitutional va- 
cancy complex, as given in the second column of Table 
IIVI show the same behavior. This makes it plausible that 
the 5sp and 6sp impurities form a split-vacancy complex 
with the vacancy, while the 4sp impurities prefer the sub- 
stitutional complex. The figure also shows the unusual 
behavior of the Al impurity, leading to a relatively large 
outward relaxation in Si and only a small inward relax- 
ation in Ge. Thus the calculations clearly indicate that 
the size of Al is larger than the one of Ga, which points to 
an unusual behavior of Al, showing also up in the stabil- 
ity of the Al-split-vacancy complex in Si-bulk. Thus the 
NN relaxations of the isolated impurities correlate qual- 
itatively well with the stability of the impurity-vacancy 
complexes. However a detailed comparison is not possi- 
ble. For instance, the outward NN relaxations of the Al 
impurity in Si are slightly smaller than the ones of As; 
yet Al forms with the vacancy in Si a split complex, but 
As a substitutional one. Similarly, the NN relaxations of 
the Si impurity in Ge are inwards, but in the substitu- 



tional Si-vacancy complex Si relaxes by 12% towards the 
vacancy. 



VI. SUMMARY 

In this paper we have studied the electronic and geo- 
metrical structure of vacancy complexes with "oversized" 
and "normal-sized" impurities in Si and Ge. First we dis- 
cuss the electronic structure of the Sn- vacancy complex. 
In agreement with recent ab initio studies for the Sn- 
vacancy complex in S:^'"* and in analogy to our recent 
studies^ for the Cd-vacancy complex in Si and Ge, we 
find a split-vacancy complex as the stable configuration 
in both hosts with the impurity on the bond center po- 
sition and the vacancy split into two "half vacancies" on 
the neighboring sites. The density of states of Sn and Cd 
are characterized by the same sequence of states and sim- 
ilar charge states exist. The calculated hyperfine fields 
of Sn and the Si neighbors are in good agreement with 
the measurements of Watkinai; however, no agreement is 
obtained with the available isomer shift data for Sn^si. 

In the second part of the paper we present a systematic 
study of vacancy complexes with 11 different impurities, 
both in Si and Ge. We find that, in both hosts, impurities 
of the 5sp and 6sp series form split-vacancy complexes, 
while impurities of the 4sp and 3sp series prefer, more 
or less distorted, substitutional complexes. An excep- 
tion from this rule is Al, forming in Si a split complex 
and in Ge a strongly distorted substitutional complex. 
Qualitatively we explain the results in terms of the size 
of these impurities, such that oversized impurities can 
lower their energy in the "free space" available at the 
bond center site. To examine the "size" of these impuri- 
ties, we calculated the nearest neighbor relaxations of the 
isolated impurities and find a good correlation between 
the calculated NN relaxations and the structure of the 
impurity- vacancy complexes. 
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